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m . ABSTRACT 

CN| ■ We show that the standard model for geometrically thin accretion disks (a-disks) 

leads to inconsistencies if selfgravity plays a role. This problem arises from the 
parametrization of viscosity in terms of local sound velocity and vertical disk scale height. 
A viscosity prescription based on turbulent flows at the critical effective Reynolds num- 
ber allows for consistent models of thin selfgravitating disks, and recovers the a-disk 
solution as the limiting case of negligible selfgravity. We suggest that such selfgravitat- 
ing disks may explain the observed spectra of protoplanetary disks and yield a natural 
explanation for the radial motions inferred from the observed metallicity gradients in 
disk galaxies. 
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1. Introduction 

One of the major shortcomings of the current the- 
oretical descriptions of accretion disks is lack of de- 
tailed knowledge about the underlying physics of vis- 
cosity in the disk. This problem is significant because 
almost all detailed modelling of the structure and evo- 
lution of accretion disks depends on the value of the 
viscosity and its dependence on the physical param- 
eters. While there is general agreement that molecu- 
lar viscosity is totally inadequate and that some kind 
of turbulent viscosity is required, there is far less 
certainty about its prescription. Most investigators 
adopt the so-called a-ansatz introduced by Shakura 
(1972) and Shakura and Sunyaev (1973) that gives 
the viscosity [v) as the product of the pressure scale 
height in the disk (h) 1 the velocity of sound (c s ), and 
a parameter a that contains all the unknown physics. 
One interprets this as some kind of isotropic tur- 
bulent viscosity v = vt = hvt where h is an (un- 
known) length scale and Vt an (unknown) character- 
istic velocity of the turbulence. One may then write 
a = (v t /c s ) ■ (h/h). On general physical grounds 
neither term in parentheses can exceed unity so that 
a < 1. If initially v t > c s , shock waves would result in 
strong damping and hence a return to a subsonic tur- 
bulent velocity. The condition l t > h would require 
anisotropic turbulence since the vertical length scales 
are limited by the disk's thickness, which is compara- 
ble to h. 

A parameterization of this sort for v is of course 
only of interest if the proportionality parameter, a, is 
(approximately) a constant. One can expect this to 
happen only if the scaling quantities are chosen in a 
physically appropriate manner. Models for the struc- 
ture and evolution of accretion disks in close binary 
systems (e.g., dwarf novae and symbiotic stars) show 
that Shakura and Sunyaev's parameterization with a 
constant a leads to results that reproduce the overall 
observed behaviour of the disks quite well. As a result 
of this success, the a-ansatz is now used in practice in 
all sorts of accretion disks. It is noteworthy, however, 
that the a-ansatz retains no information about the 
mechanism generating the turbulence but only about 
limits to its efficiency. 

In this contribution, we wish to 

• show that simple application of the above vis- 
cosity prescription to the case of selfgravitating 
accretion disks leads to an apparent internal in- 



• propose a possible resolution of this problem in 
terms of a viscosity prescription (Sec. ||). 

• discuss the implications of the new formulation 
of viscosity for the structure and evolution of 
accretion disks (Sec. Q). 

• discuss protostellar and galactic disks as two ex- 
amples where a non-standard viscosity descrip- 
tion in selfgravitating disks is important (Sec. 

D- 

2. The Structure of Selfgravitating Accretion 
Disks 

2.1. Conditions for Selfgravity in Accretion 
Disks 

In the following we assume that the accretion disks 
are geometrically thin in the vertical direction and 
symmetric in the azimuthal direction. We approx- 
imate the vertical structure by a one zone model. 
Then a disk model is specified by the central mass 
M*, the radial distributions of surface density S(s), 
central plane temperature T c (s) and effective temper- 
ature T e f{(s). The relevant material functions are the 
equation of state, the opacity and the viscosity pre- 
scription. 

One can estimate the importance of selfgravity by 
comparing the respective contributions to the local 
gravitational accelerations in the vertical and radial 
directions. 

The vertical gravitational acceleration at the disk 
surface is 2wGT, and GM^h/s 3 , respectively, for the 
selfgravitating and the purely Keplerian case, respec- 
tively. Selfgravitation is thus dominant in the vertical 
direction when 



M d 7ts 2 £ 1 h 
M* ~ M, > 2? 



(1) 



where M d (s) is the enclosed mass in the disk and is 
given approximately by M d ~ 7rs 2 S. With character- 
istic numbers for the relative geometrical thickness 
of non-selfgravitating disks of h/s ~ 0.1 ... 0.02, 
condition (|lj) translates into a condition for vertical 
selfgravity of M d /M* > 0.01 . . . 0.05. 

Similar considerations lead to the condition 



Ma > M* 



(2) 



consistency (Sec. 2.2). 



for selfgravitation to dominate in the radial direction. 
Thus, for increasing disk masses selfgravity first be- 
comes important in the vertical direction. 

2.2. Keplerian Selfgravitating Disks (KSG) 

In this Section, we describe a selfgravitating (SG) 
accretion disk in the framework of the model by 
Shakura and Sunyaev (1973; hereafter referred to as 
the standard model), changing only the equation of 
hydrostatic equilibrium in the direction perpendicu- 
lar to the disk's plane (vertical or z direction) to ac- 
count for selfgravity in the direction perpendicular to 
the disk. Thus, while in the standard model the local 
vertical pressure gradient is balanced by the z com- 
ponent of the gravitational force due to the central 
object, in the SG case we have balance between two 
local forces, namely the pressure force and the grav- 
itational force due to the disk's local mass. For the 
present purpose, we consider only selfgravity in the 
vertical direction, and assume that in the radial di- 
rection centrifugal forces are balanced by gravity from 
a central mass (Keplerian approximation). In the fol- 
lowing, we will often refer to such accretion disks, in 
which selfgravity is important only in the vertical di- 
rection, as Keplerian selfgravitating (KSG) disks. 

Hydrostatic equilibrium in the vertical direction 
yields 

P = ttGS 2 (3) 

(Paczyhski 1978), where P is the pressure in the cen- 
tral plane (z = 0), £ is the surface mass density inte- 
grated in the z direction, and G is the gravitational 
constant. Equation (H) is strictly valid only in the 
limit of dominating selfgravity. We will comment on 
the case of marginal selfgravity later. To deduce Eq. 
rtS), we have made a one- zone approximation for the 
disk's vertical structure, i.e., we have replaced deriva- 
tives with respect to the z direction by quotients of 
differences between the disk's surface and the central 
plane. 

Since details of the thermodynamics in the z direc- 
tion are of no particular relevance to our argument, 
we shall likewise assume the disk to be isothermal in 
the vertical direction. 

Integrating the equation of conservation of angular 
momentum gives 



vY> 



M 



with the radial mass flow ratey M, the Keplerian ro- 
tational frequency u), its radial derivative a/, and a 
quantity £ allowing for the integration constant or, 
equivalently, for the inner boundary condition. For 
a detailed discussion of £ see, e.g., Duschl & Tschar- 
nuter (1991), Popham & Narayan (1995) and Donea 
& Biermann (1996). For simplicity, we set the bound- 
ary condition £ = in the subsequent discussion. This 
does not alter the essence of our argument, and only 
changes details close to the disk's inner radial bound- 
ary, since the product s 2 u> increases with s. 
Finally, we use the standard viscosity ansatz 



v = ahc s 
and have for the sound velocity 



P 



(5) 



(6) 



Equations (§|) and (|) give h = c 2 /(27rG£). Substi- 
tuting into Eq. (||), we get uH = (ac\) / (2irG) . Com- 
paring this with Eq. (H), we find that the set of Eqs. 
(pf) - (0) can be fulfilled only for a certain sound ve- 
locity, or temperature given by: 



c? = - 



and hence 



2ttGM lu 2-kGM ( dlogs 

dlogw 



kT c (2ir.. 

— - = — GM 
niH V 3a 



2/3 



(7) 



(8) 



Thus, for a KSG disk, the a-ansatz leads to the re- 
quirement of a constant temperature for all radii s 
(or, if £ ^ 0, the temperature is prescribed as a func- 
tion of s), independent of thermodynamics. It also 
requires that the disk structure satisfy 



h£ = 



M 2 ' 3 



2ttG 



(9) 



(3a) 2/3 (27rG) 1/3 ' 

In a standard non-SG accretion disk the temperature 
is a free parameter which is determined by the energy 
released by the inward flow of the disk gas (M), by the 
local viscosity, and by the respective relevant cooling 
mechanisms. The viscosity depends on T c via Eq. 



7 (^-0 



(-1) 



2 We choose the convention of radial mass flow rate M and radial 
velocity v s positive for inward motion. 



(||) and on the equation of hydrostatic support in the 
direction normal to the disk, namely 



v,p 



(10) 



(v$ = suj) which in the non-SG case replaces Eq. (Q). 

In the KSG case, it is the surface density £ (and 
hence h) which must adjust in order to radiate the en- 
ergy deposited by viscous dissipation and provided by 
the inward flowing material. While detailed solutions 
are beyond the scope of this paper, they clearly exist 
formally. On the other hand, the normal thermostat 
mechanism does not operate, at least in the steady 
state. Indeed in certain circumstances, the condition 
of constant mid-plane temperature appears to be in- 
consistent with the basic thermodynamic requirement 
that the average gas temperature in the disk exceed 
that of the black body temperature required to ra- 
diate away the energy dissipated by viscous stresses 
(see Appendix). We are therefore doubtful whether 
a physically plausible and stable quasi-steady state 
solution exists. 

Furthermore, the Jeans condition for fragmenta- 
tion in the disk into condensations of radius R is 



— qOpR 2 > c 2 



(11) 



(see Mestel 1965) where q is factor of order unity. On 
the other hand, Eqs. (||) and (^) give 



2kGY.1i = 4nGphr = a 



(12) 



Thus, a selfgravitating disk is on the verge of frag- 
menting into condensations of radius R ~ h unless 
these are destroyed by shear motion associated with 
the Keplerian velocity field. 

Finally we note that the Reynolds number K = 
vl/fmoi in the disk flow is extremely high in any as- 
trophysical context and this in itself is likely to lead 
to strong turbulence. 

Indeed, as has been pointed out for example by 
Lynden-Bell and Pringle (1974) and by Thompson et 
al. (1977), the high Reynolds number character of 
the flow will lead to the generation of turbulence and 
hence to a steady enhancement in the effective (or 
Eddy) viscosity until the Reynolds number becomes 
subcritical. The limiting viscosity in this case is given 
for a disk as 



SRcrit 

^10 2 



Vrj) s = (3v (jl s (13) 

. . . 10 3 . This form is quite 



where 5ft crit ~ 1//3 
different to the a-ansatz. 

For all of the above reasons, we are led to consider 
alternative prescriptions for viscosity which may be 
more relevant to selfgravitating Keplerian (or fully 
selfgravitating) disks. 

3. Prescription for Turbulent Viscosity 



As noted in Sections [j] and 2.2 the need for some 
kind of turbulent viscosity in accretion disks is gener- 
ally recognized, as is the very high Reynolds number 
of the flow in the absence of such a viscosity. In any 
event, it seems reasonable to assume that the turbu- 
lence is driven by the velocity field in the disk which 
has characteristic length and velocity scales s and ty,, 
respectively. We might then expect that 



Vth oc Av^As 



(14) 



where Av,p and As are representative velocity and 
length scales of the flow. Furthermore, we may write 



Ai^ 



dv 6 «0 

-As ~ —As 



(15) 



ds s 

If we consider turbulent elements in a smoothed out 
background gas with sound speed c s , we may impose 
the limit that the turbulent velocities do not exceed 
c s . Thus Eq. (M) gives 



—As ~ Atu < c s 
s 



As 



s s 

— Av^ < — c s 

v<p V<p 



(16) 



(17) 



From Eq. (10), this implies that for a standard Kep- 
lerian disk, 



and hence that 



As < h 



ihc s 



(18) 



(19) 



where a < 1. Note that the upper bound to As im- 
plies approximately isotropic turbulence. This is the 
standard a-ansatz but derived from considerations 



of rotationally generated turbulence. The argument 
can, in fact, be inverted to show that under assump- 
tion of isotropic turbulence Eqs. ( JIOJ ) and ([19), with 
As = h, imply Av^ ~ c s which again yields the func- 
tional dependence of the a-ansatz. Implicit in this 
argument is that the effective viscosity is determined 
by the largest physically reasonable length and veloc- 
ity scales. 

For a SG disk (whether Keplerian or not) Eq. ( |l0| ) 
is no longer valid so this analysis does not apply. In 
physical terms, the scale height in the disk no longer 
reflects global properties of the disk (mass of and dis- 
tance to the central star) but is set by local condi- 
tions. Furthermore, if the material is significantly 
clumped, as seems likely, for example, in molecular 
clouds in galactic disks, an alternative prescription 
may be more appropriate. We suggest that in such 
situations (and possibly in others also) the appropri- 
ate length and velocity parameters may be the radial 
distance, s, and the azimuthal velocity, v$ or, more 
specifically, some fractions of each, since in some sense 
they already represent the maximum of either scale 
available. 

In support of this choice of s as the natural length 
scale we note that (a) it is the only length scale which 
is relevant for angular momentum transport which 
contains information about the driving agent for the 
turbulence - namely the rotation field, and (b) it is 
relevant to several other processes for transporting 
angular momentum such as spiral arms and large scale 
magnetic fields. Likewise, the orbital velocity is the 
only relevant velocity scale containing information on 
the rotation. 

In these circumstances, we should obtain the Rey- 
nolds form for the viscosity, namely 



v = fR = Psv^ 



(20) 



where (3 is a constant. Since the rotation field is as- 
sumed to be driving the turbulence, it seems appro- 
priate to restrict (3 using the condition that the effec- 
tive Reynolds number of the flow not fall below the 
critical value 5ft cr ;t for the onset of turbulence. Thus 
(3 should satisfy 



0< 



1 



% : 



10" 



10" 



(21) 



In the following we will refer to Eq. (E0J) as the 0- 
ansatz and to the disk structure arising from this 
viscosity prescription as (3- disks. We suggest that 



this formulation is most appropriate for selfgravitat- 
ing disks which we explore further in Sec. 0. 

We note, however, that if the disk matter distri- 
bution is clumpy (e.g., clouds within a low density 
smoothed out distribution) then there is a formal con- 
nection between the a- and the /3-prescriptions. Since 
in the /3-formulation the clump velocities are of order 
«0 shock heating will tend to heat the low density 
interclump gas until its sound speed c s ~ v^. The in- 
terclump gas will then have a scale height h ~ s, the 
scale of the clumpy disk, and will hence be roughly 
a spherical structure. At this point the a- and (3- 
prescriptions look formally identical but the scale 
height and sound speed now refer to a more or less 
spherical background distribution of hot gas in which 
a disk structure of cloudy clumps is inbedded. 

4. Structure of Selfgravitating f3 Disks 

For SG /3-disks, Eqs. (g), (@) and (§) are applicable, 
but Eq. (JsJ) must be replaced by Eq. (pfj|). It then 
follows that 



M 



and that 



nGh M 



irGhM 



(22) 



(23) 



v slo' (3s 3 lu' 

Thus the SG /3-disks recover the thermostat property 
of the standard disk, namely that the temperature 
and scale height can adjust to accomodate (radiate 
away) the energy input to the system from viscous 
dissipation and inward motion. 

4.1. Keplerian Selfgravitating Disks 

For the particular case of a KSG /3-disk we have: 



h 



4irGM 
3" 



f 



(24) 



(3 (GM,)!/2 s l/2 

For the SG /3-disk it follows immediately from Eq. 
(|22j) and from mass conservation in the disk that the 
radial inflow velocity v s is given by 



M f3s 2 uj' 

Vs ~ 2tt.sE ~ 2^~ 
For the KSG /3-disk, we then have 



(25) 



-J3suj = —fivcj, 



(26) 



47T 47T 

Thus at each radius the inward velocity is the same 
fraction of the local orbital velocity. From Eq. ( |25| ) 
this, in fact, holds for any SG /3-disk in which the 
angular velocity is a power law function of s with 
adjustment only to the numerical factor in Eq. ([26|). If 
satisfies the constraint (£lj) , then the approximation 
of centrifugal balance in the radial direction remains 
well justified. 

Under these conditions the dissipation per unit 
area of a SG /3-disk is given by 



D 




2^ c 4 ff 



(27) 



where a is the Stefan-Boltzmann constant. For an 
optically thick KSG disk this yields the same radial 
dependence of T e g as for the standard disk, namely 



Tf.fr = 



( GMM* 



1/4 



-3/4 






87TO- 



(28) 



This temperature dependence which is identical to 
that of the standard model then leads to the well 
known energy distribution for an optically thick stan- 
dard disk of S v oc i/ 1 / 3 . This also implies that - as 
long as the disks are not fully self-gravitating - it is 
hard to distinguish between an a- and a /3-disk model 
observationally. 

4.2. Fully Selfgravitating Disks (FSG) 

We turn now to the case of the fully selfgravitating 
(FSG) /3-disk, in which the disk mass is sufficiently 
great that it dominates the gravitational terms in 
the hydrostatic support equation in both the radial 
and vertical directions. While there are many po- 
tential solutions for the FSG disk structure, one is 
well known in both mathematical and observational 
terms, namely the constant velocity (w^ = const.) 
disk. Within such a disk structure we have simulta- 
neous solutions to the equation of radial hydrostatic 
equilibrium and Poisson's equation of the form 



«0 = soj = vq and SocS 



•so 



(29) 



(Toomre 1963; Mestel 1963). For the FSG disk, Eq. 
(22) then leads to 



M 
f3sv 



(30) 



which has the same radial dependence as the struc- 
tural solution as shown in Eq. (p9|). Thus Eq. (30) 
may be viewed as giving the rate of mass flow through 
the disk for a FSG /3-disk with constant rotational ve- 
locity vq. Finally, the equation of continuity provides 
a constraint if the structure is to maintain a basically 
steady state structure. We then have: 



5E d, 3d 



0, Os—toBs&W (31) 

For the v^ = vq = const, disk, it then follows from 
Eq. @ that dY./dt = 0. Thus the constant velocity 
disk represents a steady state solution in regions suf- 
ficiently far from the inner and outer boundaries of 
the /3-disk. 

It is then possible, in the spirit of the discussion of 
Eqs. (|27| ) and (|2g|), to calculate the energy dissipation 
rate per unit area D for the constant velocity /3-disk. 
We then find 



D = 2aT c % = 



Mvl 
Airs 2 



so that 



T rt 



Mvl 



1/4 



-1/2 



(32) 



(33) 



For flux density S v emitted by an optically thick, con- 
stant velocity /3-disk it then follows that 



S v oc v 



(34) 



In reality, a sufficiently massive disk may be expected 
to have an inner Keplerian (standard) zone, a Ke- 
plerian selfgravitating zone (KSG), and a fully self- 
gravitating zone (FSG). We should therefore expect 
a smooth transition in the spectral energy distribu- 
tion from the v 1 ' 3 spectrum of the inner two zones to 
the v~ l spectrum arising at longer wavelengths from 
the FSG zone. One could turn this argument around 
and argue that, if no other components contribute to 
the spectrum, the flatness of the vF v distribution is a 
measure for the importance of selfgravity and thus for 
the relative mass of the accretion disk as compared to 
the central accreting object. This, of course, applies 
only to the optically thick case which may not arise 
frequently in strongly clumped disks. 



4.3. Time Scales 

The evolution of accretion disks can be described 
by a set of time scales. For our purposes, the dy- 
namical and the viscous time scale are of particular 
interest. 

The dynamical time scale Td yn is given by 



1 

7"dyn — • 



(35) 



While this formulation applies to all cases, selfgrav- 
itating or not, it is only in the non-SG and in the 
KSG cases that ui is given by the mass of the central 
accretor and by the radius. In the FSG case, u> is 
determined by solving Poisson's equation. 

The timescale of viscous evolution r v i sc is given by 



Tvisc = — (36) 

v 

In the standard non-SG and geometrically thin (h <C 
s) case, this leads to 



non— SG 



7"dyn ^^ 7"dyn 



( SV ^oyn > ^ (37) 

\hJ a a 



In KSG and FSG disks, r visc is given by 



KSG 



FSG 



_SG 



7*dyn 



(38) 



With a < 1 and /3 < 1 (Eq. gTj) under all cir- 
cumstances T V i SC 3> Tdyn- In the SG cases the ratio 
between the two timescales decouples from the disk 
structure. In all cases the models are self consistent in 
assuming basic hydrostatic equilibrium in the vertical 
direction. 

5. Possible Applications 

5.1. Protoplanetary Accretion Disks 

T Tauri stars have infrared spectral energy dis- 
tributions vF v which can be approximated in many 
cases by power laws vF v oc v n with a spectral index 
n in the range ~ ... 1.3. Assuming this spectral be- 
haviour to be due to radiation from an optically thick 
disk, it translates into a radial temperature distribu- 
tion T c g oc s~ q with n = 4 — 2/q. 

A passive disk that only reradiates reprocessed ra- 
diation from the central star has q = 3/4 or n = 4/3. 
An optically thick active non-selfgravtating accretion 
disk which radiates energy that is liberated through 



viscous dissipation shows the same spectral distribu- 
tion. Thus, if accretion disks were the only contrib- 
utors to T Tauri spectra, and if the disks were opti- 
cally thick and non-selfgravitating, the spectral shape 
should be the same for all objects. This is in clear con- 
trast to the above mentioned observed spectral energy 
distributions. 

Adams, Lada and Shu (1988) were the first to dis- 
cuss the possibility of a non-standard radial temper- 
ature distribution with q =/= 3/4. Using q as a free 
parameter, they find that for flat spectrum sources, 
their best fits require disk masses that are no longer 
very small as compared to the accreting stars' masses. 
They already mention the possibility that the flatness 
of the spectrum and selfgravity of the disk may be re- 
lated. On the other hand, at that time this indirect 
argument was the only evidence for large disk masses. 

Due to the lack of accretion disk models yielding 
the appropriate values of q, since the work of Adams, 
Lada and Shu several other interpretations for the 
flatness of the spectra were offered, for instance dusty 
envelopes engulfing a star with a standard disk around 
it (Natta 1993). 

In the meantime, high resolution direct observa- 
tions of protostellar disks yield independent strong 
evidence for comparatively large disk masses. Lay et 
al. (1994), for instance, find a lower limit for the disk 
masses in HL Tau - one of the sources in Adams, Lada 
and Shu sample of flat spectrum T Tauri stars - of 
-O.O2M . 

We suggest that the flatness of the spectrum actu- 
ally reflects the mass of the disk, i.e., the importance 
of selfgravity. For disk masses considerably smaller 
than ~ 1/3QM*, the standard accretion disk mod- 
els apply. For disks whose masses are larger but still 
small compared to M* the spectral behaviour is not 
altered significantly, but disk structure and the time 
scale of disk evolution (t v ; sc , see Eqs. [37] and [38]) 
change. For even more massive disks, we expect a 
clear trend towards flatter spectra that approach an 
almost constant vF v distribution if selfgravity in the 
disk becomes important in radial as well as in verti- 
cal direction. But the time scale for the disk evolution 
remains the same (Eq. |38| ). 

5.2. Galactic Disks 

The relevance of viscosity in the evolution of galac- 
tic disks has been the subject of discussion since von 
Weizsacker (1943, 1951) and Lust (1952) first raised 



the issue nearly fifty years ago. They noted then that, 
with an eddy viscosity formulation (a /3-disk), the 
time scale for evolution of typical galactic disks was 
comparable to the age of the universe and suggested 
that this might account for the difference between spi- 
ral and elliptical galaxies. 

With the subsequent realization that galactic disks 
moved primarily under the influence of extended mas- 
sive halos, interest in FSG disks waned. However, as 
noted above, it is possible for a massive disk to exist 
and evolve under the influence of viscosity embedded 
in such a halo gravitational field. Indeed, in the event 
that such a structure forms, it must evolve under vis- 
cous dissipation and can achieve a quasi-steady state 
with essentially the same mass and energy dissipation 
distribution as for the FSG constant velocity disk. We 
refer to this case as an Embedded Self Gravitating 
(ESG) disk. 

The timescale for viscous evolution t v - 1sc as given in 



section 4.3 suggests a means of differentiating between 
the a- and /3-formulations for this case. For a normal 
spiral galaxy with a suggested mean temperature in 
the gaseous disk of around 10 4 K and a scale height 
of around 300 pc, we obtain 



r («) 

vise 

M 



10*Tdyn~310 u yr 



KT 



Krr dyn - 3 10 y - 3 10 1U yr (39) 



Thus, with these parameters, little evolution would 
take place in a Hubble time on the a-hypothesis but 
significant evolution is predicted on the /3-hypothesis. 

In terms of inflow velocities the /3-ansatz suggests 
values in the range 0.3 — 3kms _1 which would be 
exceedingly hard to measure directly: the a-ansatz 
suggests still lower values. On the other hand, it may 
be possible to provide limits on the viscosity through 
other observational constraints. For example, the 
build up of the 3 kpc molecular ring in our own galaxy 
can be interpreted as due to viscosity driven inflow in 
the constant velocity part of the galactic disk which 
ceases (or at least slows down) in the constant an- 
gular velocity inner regions (Icke 1979; Dather and 
Biermann 1990). Similarly, several authors have sug- 
gested that the radial abundance gradients observed 
in our own and other disk galaxies may be due to 
radial motion and diffusive mixing associated with 
the turbulence generating the eddy viscosity (Lacey & 
Fall 1985; Sommer-Larsen and Yoshii 1990; Koppen 
1994; Edmunds & Greenhow 1995, Tsujimoto, Yoshii, 



Nomoto and Shigeyama 1995). According to these au- 
thors, radial inflows of around lkms -1 at the galac- 
tic location of the Sun are required for optimum fits 
to the abundance gradient data within the context 
of the viscous disk hypothesis. Such inflow velocities 
are consistent with the /3-ansatz but could, of course, 
be generated also by other means (e.g effects of bars, 
magnetic fields). 

6. Summary 

The standard model for geometrically thin accre- 
tion disks with viscosity proportional to sound ve- 
locity and vertical scale height (often referred to as 
a-disks) leads to inconsistencies if the disk's mass 
is large enough for selfgravity to play a role. This 
problem arises even in Keplerian selfgravitating disks 
in which only the vertical structure is dominated by 
selfgravity while the azimuthal motion remains Kep- 
lerian. 

We propose a viscosity prescription based on the 
assumption that the effective Reynolds number of the 
turbulence does not fall below the critical Reynolds 
number. In this parametrization the viscosity is pro- 
portional to the azimuthal velocity and the radius (/?- 
disks). This prescription yields physically consistent 
models of both Keplerian and fully selfgravitating ac- 
cretion disks. Moreover, for the case of sufficiently 
small disk masses, we recover the a-disk solution as a 
limiting case. 

Such /3-disk models may be relevant to protoplan- 
etary accretion disks as well as for galactic disks. In 
the case of protoplanetary disks they yield spectra 
that are considerably flatter than those due to non- 
selfgravitating disks, in better agreement with ob- 
served spectra of these objects. In galactic disks, they 
result in viscous evolution on time scales shorter than 
the Hubble time and thus offer a natural explanation 
for an inward flow that could account for the observed 
chemical abundance gradients. 
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A. Thermodynamic Considerations for KSG 

a-disks 

For a KSG a-disk, we have from Eq. H that 



(Al) 
where m is the mass flow rate in solar masses per year, 
and a_i = a/0.1. 

If the disk is optically thick and advection is negli- 
gible, viscous dissipation leads to local effective tem- 
perature of 



T c e — 



{GM) l/A MV*8-*'* 

1/4 



.5310 J K 



(A2) 



with to the mass of the central star in solar units and 
sa the radius in astronomical units. 

An essential thermodynamics requirement is that 
T c > T c g or that 
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— = 3.53 10 — A a\ < 1 (A3) 

c S A 



This condition is satisfied provided that 

3/5 8/5 

A m ' a 1 
to > tot = 3.27 10" 
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(A4) 



and that the disk is selfgravitating in the vertical di- 
rection at sa- The latter condition leads to a second 
requirement on m. 

For a standard Keplerian disk, the mass flow rate 
is given (Eqs. Ifippj) h ? 



M d M A v fh\ 

M « — - w — %- = aM d ( - ) lu (A5) 



with Ma the disk's mass. From Eq. El the condition 
that the disk is non-selfgravitating is M d < {h/2s)M Sf 



and hence, from Eq. A5, that 



M < 



n fh\ 3 fGM^ 1 ' 2 

\ s 3 



(A6) 



m < to g = 3.1410 x a_i 



3 3/2 
TO J/ 



3/2 



(A7) 



A selfconsistent and physically acceptable solution 
can be obtained only if toq > TOt> that is the disk 
becomes selfgravitating at values of m which are suf- 
ficiently high that thermodynamic requirements are 
not violated. This condition may then be written as 



- > 1.0110 _1 m 



-1 -3/10 -1/10 1/5 
m i r ' rv ' 



(A8) 



Thus thin KSG a-disks with to ^> 1, sa ^> 1 appear 
to be inconsistent with basic thermodynamic require- 
ments. 

No such inconsistencies occur for /?-disks for which 
the normal thermostat effect is free to operate. 
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